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Abstract 

We construct the explicit Q-operator incorporated with the s?2-loop-algebra symmetry of 
the six-vertex model at roots of unity. The functional relations involving the Q-operator, the 
six-vertex transfer matrix and fusion matrices are derived from the Bethe equation, parallel to 
the Onsager-algebra-symmetry discussion in the superintegrable iV-state chiral Potts model. We 
show that the whole set of functional equations is valid for the Q-operator. Direct calculations 
in certain cases are also given here for clearer illustration about the nature of the Q-operator in 
the symmetry study of root-of- unity six- vertex model from the functional-relation aspect. 



1999 PACS: 05.50.+q, 02.20.Tw, 75.10Jm 

2000 MSG 17B65, 39B72, 82B23 

Key words: Six-vertex model, s^-loop algebra, Chiral Potts model, Bethe ansatz 



1 Introduction 

The aim of this paper is to construct the Q-matrix for the six-vertex model at roots of unity 
incorporable to the s/2-loop-algebra symmetry found in [14], and the functional relations proposed 



1 



in [19, 29]. Since initiated by Bethe in 1931 [13], the free energy of the six- vertex model and the 
spectrum of the XXZ-spin chain with the periodic boundary condition, 

-1 L 

where a\ is the Pauli spin matrix at the site i with L + 1 = 1, have long been studied by means 
of Bethe ansatz. The results obtained from Bethe ansatz were shown by Baxter in the seventies to 
follow from the TQ-relation method by introducing an auxiliary 'Q'-operator, a genuine invention 
in his solution of the eight-vertex model [2, 3, 4, 5]. The TQ-relation method enables one to study 
the eigenvalues of the transfer matrix T without knowing the eigenvectors, but instead to work on 
the eigenvalues of the Q-operator. In this manner, there exist many Q-matrices which satisfy the 
TQ-equation (see, e.g. [5, 23] and references therein), but only a certain special one is expected to 
exhibit the additional symmetry appeared in the six-vertex model when the anisotropic parameter 
q with A = 7j(q + c/ -1 ) is a root of unity. In this work, we find a such Q-operator. We demonstrate 
first by direct calculations in examples, then by means of mathematical arguments in general cases, 
that the Q-operator we have constructed here is the one for the symmetry study of the root-of-unity 
six-vertex model, much in the same way as the chiral Potts transfer matrix in the Onsager-algebra 
symmetry study of the superintegrable iV-state chiral Potts model (CPM) in [28]. 

In their study of the root-of-unity symmetry in eight- vertex model [19], Fabricius and McCoy 
proposed the conjectural functional relations for the eight-vertex transfer matrix T and Q72-matrix 
in [3], analogous to the set of functional equations known in the A^-state CPM [11]. After the 
occurrence of r^-de generacy in CPM was found to appear only in the superintegrable case [27], 
the Fabricius-McCoy comparison between the root-of-unity eight-vertex model and superintegrable 
CPM was further analysed, and their common mathematical structure led to the discovery of the 
Onsager-algebra symmetry of r^-model in the superintegrable CPM [28]. These exact results in 
CPM can serve as a valuable scheme in the study of symmetry problems of solvable lattice models, 
among which the root-of-unity six-vertex model is a distinguished one, due to the importance 
of the related XXZ-spin chain (1.1), and the progress made on the newly found s/2-loop- algebra 
symmetry in roots of unity case [14, 15, 16, 17, 18]. Indeed, the Fabricius-McCoy comparison in 
the limiting case of eight-vertex model with the vanishing elliptic nome strongly suggested that 
a unified theory exists for both the root-of-unity symmetry of six-vertex model and the Onsager- 
algebra symmetry of CPM from the functional relation aspect, which was clearly spelled out in 
[29]. Among the functional relations, the QQ-relation, equivalent to the Q-functional equation, 
encodes the essential features reflecting the symmetry of the model. In the theory of CPM, the 
QQ-relation (called the TT-relation in [11]) played the vital role in the derivation of the whole 
set of functional equations, by which the solution of eigenvalue spectrum of the CPM transfer 
matrix, considered as the Q-operator, was obtained [6, 25]. Works in [1, 7, 8] about the Q- 
eigenvalues in the superintegrable CPM has paved the way to both the qualitative and quantitative 
understanding of the degeneracy and Onsager-algebra symmetry of r^-models [28]. These results 
encourage us to view the CPM as a useful 'toy'-model to illustrate the symmetry nature of the 
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root-of-unity six-vertex model, and hopefully to other models. For the purpose of displaying the 
common symmetry characters governing various models we will consider in this scheme, handling 
the precise formulation which correctly and uniformly presents the functional relations for all those 
models is a non-trivial problem. By adjusting some suitable normalizing factors to the original 
operators appeared in the CPM functional relations in [11], we refine the quantity to extract a 
(somewhat unpleasant or more complicated) representation of functional relations, which can be 
extended successfully to include both the root-of-unity six- and eight-vertex models so that the 
symmetry property of the Q-operator will be clearly presented. Furthermore, we also hope the 
effort of understanding the common integrable structure of those models will help to provide some 
clues to solutions of certain unsolved important problems in CPM, such as the CPM correlation 
functions, parallel to the results known for the correlation functions of the six- and eight- vertex 
models. In the context of the -model where the six- vertex model is with a particular field (see 
[10], page 3), Bazhanov and Stroganov showed that the column transfer matrix of the L-operator 
((2.19) of [12])) possesses the properties of Baxter's Q-matrix with features of symmetry operators 
similar to, but not exactly the same as, the s/2-loop algebra which plays the role of symmetry 
in the 'zero-field' six- vertex model ([9] Sect. 3). Indeed, the Onsager-algebra symmetry operators 
inherited from the Baxter's Q-matrix provide the precise description about the symmetry structure 
of T^-model in the superintegrable CPM [28]. In the study of 'zero-field' six- vertex model, which 
we loosely call the six-vertex model in this paper, the explicit form of the six-vertex Q-operator 
which displays the root-of-unity symmetry is our main concern. This Q-operator will behave like 
the (^-operator in the root-of-unity eight- vertex discussion in [19], which was invented by Baxter 
in 1972 for his study of eigenvalues of the eight- vertex transfer matrix [3]. In this work, we produce 
a Q-operator of the six- vertex model at roots of unity by following Baxter's method of constructing 
Q72 in [3], and demonstrate that this special Q-operator is indeed in accordance with the root-of- 
unity symmetry of six- vertex model from the functional-relation aspect. Note that in the process 
of constructing this six-vertex Q-operator, we first produce certain operators Qr,Ql, then Q by 
a normalized factor as in the eight- vertex case [3]. The Qr,Ql here cannot be reached from the 
eight-vertex Qr- and Qi-operator by taking the vanishing elliptic nome limit. However in cases 
when Q72 is explicitly known, the limit exists. It should be interesting to study the comparison of 
the Q-operator in this paper with a six-vertex limit of the eight-vertex Q72 because the conclusion 
about the Q-operator in other cases not discussed in this work depends on it. Indeed the Q-operator 
we obtain in this paper coincides with Baxter's classical result on the six-vertex limit (for a generic 
q and S z = 0) of the eight-vertex Q-operator in [4] (formula (101)). 

This paper is organized as follows. In Sect. 2, we summarize the main features of the transfer 
matrix T(z) and s/2-loop-algebra symmetry of the six-vertex model at roots of unity for the use 
of later discussions. In Sect. 3, we review the functional relations in the study of Onsager-algebra 
symmetry of superintegrable CPM and the loop-algebra symmetry of the root-of-unity six-vertex 
model in [28, 29]. Through the identical mathematical structure about the Bethe equation of these 
two models related to the evaluation polynomial in describing the -degeneracy, we derive the 
set of functional equations, based on the known scheme in the superintegrable CPM case [28], 
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and then extended to the root-of-unity six-vertex model. The TQ-, fusion, and T^Q-relation are 

the established relations, discussed in subsection 3.1; while the QQ- and Q-functional equations 

are the symmetry constraints of the models, examined in subsection 3.2. Indeed, we give a more 

elaborate discussion there about the QQ- and Q-functional equation of the root-of-unity six-vertex 

model for later uses. Sect. 4 contains the main result of this paper about the Q-operator of the 

root-of-unity six- vertex model. First in subsection 4.1, by imitating Baxter's method in [3] of 

producing (^-operator in the eight-vertex model at roots of unity, we construct a Q-operator, 

depending the variable Z2 for the TQ-relation, and commuting with the spin inversion operator 

and spatial translation. Then we make direct, but non-trivial, checks for the cases of size L < 6 in 

subsection 4.2 to verify the Q-functional relation by using the explicit form of the Q-operator. The 

computational evidence has revealed the characteristics of the Q-operator, and further enhanced 

its correct role in the symmetry study of the six-vertex model at roots of unity. In subsection 4.3, 

we show by rigorous mathematical arguments that the functional relations for a general L are valid 

for the Q-operator obtained in this paper. We close in Sect. 5 with some concluding remarks. 

Notations. In this paper, Z, C will denote the ring of integers, complex numbers respectively, 

L 

Zjv = Z/NZ and i = v^l. For a positive integers n, L, we denote by <8> C n the tensor product of 

L-copies of the vector space C n . 



2 Six-vertex Model at Roots of Unity and Loop-algebra Symme- 
try 

We start with some basic notions in the six-vertex model at roots of unity considered in the later 
sections. This summary will be rather sketchy, but also serve to establish the notation: more 
detailed information can be found in any of the standard references listed in the bibliography, such 
as [14, 15, 18] and references therein. 
It is known that the L-operator 1 
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Using the L-operator, one defines the monodromy matrix, 0^ =1 L^(z2) with L^{z2) = L(z2) at 

Here we use the convention in [18], Eq. (1.3): a — isinhi(u — i7), b — — isinh|(u + i7), c = — isinhi7, with 
the variables z — —e v ,q = — e' 7 . Note that the q here differs from that in [18] by minus sign for its connection with 
U q (sl2) as Eq. (2.3) in [15], where — e _z is equal to Z2 here. 
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the site £, and its trace 

L 

T(z) = tr avx (® L e (z*)), for z G C, (2.2) 

L „ 

form a commuting family of (® C )-operators, which define the transfer matrix (of size L with 
the periodic condition) of the six-vertex model. The XXZ Hamiltonian (1.1) is the logarithmic 
z ^-derivative of T(z) at z = where A = ^{q + q -1 )- The transfer matrix T(z) commutes with 
the z-component of total spin S z , the spin-inversion operator R, 

S* = \j^ol R = ®of, 

and the spatial translation operator Sr, which takes the jth. column to (j + l)th one for 1 < j < L 
with the identification L + 1 = 1. In the basis ®| ± l)g for the diagonal cr| where a z \ ± 1) = ±| ± 1), 
the number of down spins n is known to relate to S z by n = ^ — S z , and (— l) ra is the quantum 
number for the operator S = <S>gcrf, which commutes with T(z). 

When q is a iVth root of unity, degenerate multiplets for the spectrum of XXZ-eigenvalues 
occur with the s^-loop-algebra symmetry for the eigenspaces [14]. In later sections, we are going 
to examine this loop-algebra symmetry of the transfer matrix T(z) through the Q-operator and 
related functional relations, as in the symmetry study of superintegrable CPM. For definiteness, in 
this paper, we will confine the discussion of the six-vertex model only to the cases for even L (the 
chain-size), and odd N (the order of q). We also assume all q, uj{:= q 2 ), q? are primitive iVth roots 
of unity. Except section 3, the discussion of the transfer matrix T(z) and all relevant operators will 
be confined only to sectors S z = (mod N), unless stated otherwise, i.e., they are regarded as 
operators of V where 

V:= Yl ft = ±l- (2-3) 

J2 e Pe=0 (mod AO 

Then the s/2-loop-algebra symmetry is generated by and the operators 

S ±(N) = J2i< n <-< JN <L 1 ® ■ ■ ■ ® 1 ® a£ ® q~ aZ ® ■ ■ ■ ® q~ aZ 
®<rj ® q~ 2aZ ■ ■ ■ ® q~( N -V aZ ® af N ® 1 ■ ■ ■ ® 1, 

r ± ( JV )= J2i< n< ... <jN < L 1 ® • • • ® 1 ® <y% ® ® • • • ® q aZ 

mf 2 ® q 2aZ ■ ■ ■ ® g^" 1 )' 7 " (g) cr^ ® 1 ■ ■ ■ ® 1, 
with the relation T ± ( Ar ) = RS^R' 1 [14, 15]. 

3 Functional Relations of Superintegrable CPM and Root-of-unity 
Six-vertex Model 

Here we give a summary of functional relations in [29], which govern both cases of the superinte- 
grable CPM and the six-vertex model at roots of unity, and a more detailed account of QQ- and 
Q-functional relation in the six-vertex model case will be discussed in subsection 3.2 for later uses. 
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In the discussion of this section, the spin S z of the six- vertex model can be arbitrary, i.e., no the 
iV-multiple restriction for S z ; while in the CPM case, the order N is not required to be odd, nor 
the chain-length L being even, and a; is a iVth root of unity. 

We begin our consideration by reviewing the Onsager-algebra symmetry of superintegrable 
CPM in [28, 29]. The rapidities of iV-state CMP are represented by coordinates of the genus- 
(iV 3 - 2N 2 + 1) curve 



W : 



ka N + k'c N = d N , i kx N = 1 — k'fi~ N , 

kb N + k'd N = c N ^ { ky N = 1 - k'fi N 



where [a, b,c, d] £ P 3 , (x,y,fi) = (§, \, G C 3 , and k' is the parameter with k 2 = 1 - k' 2 / 0, 1. 



(3.2) 



Then the variables t := xy, A := fi N define the hyper elliptic curve of genus N — 1: 

Consider the following symmetries of rapidities of order N and 2: 

U : [a, b, c, d] t— ► [wo, 6, c, d], (x, y, /it) i— ► (wi, y, /i), (i, A) i— ► (wi, A); 

C : [a,b,c,d] i-> [b,a,d,c], (x, y, //) h-> (y, x, /i" 1 ), (i, A) i-> (i, A -1 ). 

The t^ 2 ) -model is constructed from Yang-Baxter solutions G p (t) for p = [a, b, c, d] £ W: 

b 2 G p it)=( b2 ~ td2X , ^ dX J Z ), teC , 
PK ' \ -t(bc-adX)Z- 1 -tc 2 +uja 2 X J 

associated to an inhomogeneous .R-matrix of the six- vertex model [12, 27], where X, Z are the Weyl 
operator of C N : X\n) = \n + 1), Z\n) = uj n \n) for n £ Zjy. The Tp^-matrix is defined by 

L 

T^{t) := tw((g) G p/ (cot)), G p/ (t) = G p (t) at site t 
l=i 

The 'Q'-operator associated to r^-matrix is the CPM transfer matrix [11, 12], which is the (<g 
C Ar )-operator defined by 

a' a' L — 

T cp (p; s)ai:::y L = n w pA^i - °i)w Pt8 (<ri - , o h o\ e z N , 

1=1 

where p, s £ 2XJ, and W PjS , W p>s are Boltzmann weights defined by the iV-cyclic vectors: ^/'^q) = 
™ d P bs-a P c s ^ W^A^L = m ^vds-d P a a ^ ^ . f h we H-known star-triangle relation. 
Then for a fixed p £ 25} , {T cp (p; s)}<, g 22J f° rm a commuting family of operators, commuting with 

the spin-shift operator X(:= YlgXi) of ((g) C^) and the spatial translation Sr. In the theory of 
CPM, there is a collection of functional relations among the CPM-transfer-matrix T cp (p;*), the 
Tp 2 ^- and fusion t p ^ -matrices, involving the rapidity-symmetries U,C in (3.2). In this paper, the 
discussion of CPM will be confined only to the superintegrable case, i.e., the rapidity p given by 
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Hp = 1, Xp = y p = r/2 with rj := (j^p-) 17 , where simplification occurs for the functional relations. 
For simplicity, the label p will be omitted for operators at the superintegrable point p, and simply 
write r( 2 )(*), T cp (*), etc. It is convenient to use the normalized CPM transfer matrix 2 : 

Q(s) = T cp (s)(l^x^ (=; e -iPg (s)); (3 3) 

iV L (l - f] 2 X s ) L (f] 2 X s ) Pa (l] 2 y s ) P bjJL s M 

where P a , P^, P M are quantum numbers in the CPM theory and P is the total moment. Introduce 
the variable t = r]~ 1 t s and write r^-matrix by T^ 2 \t) = r^ 2 ^(t s ). By explicit results about 
eigenvalues of CPM transfer matrix in [1, 7, 8], the Onsager-algebra symmetry of r^-model was 
verified rigorously in [28] through the functional relations. Furthermore, if we normalize the t^- 
and r^-matrices by 

{t) - (i-u,-^(i-^ ' {t) ~ nfc 2 -i(i-^*) £ ' 

the TQ-, and fusion relations of CPM possess the identical structure as those in the root-of-unity 
six-vertex model through the evaluation polynomial of the respective symmetry algebra. By these 
common structures, as a theory parallel to CPM, we propose in [29] the functional relations for the 
Q-operator of the six-vertex model to display the root-of-unity symmetry, now described below in 
a uniform manner. 

3.1 Bethe equation, TQ-relation, and fusion relation 

With the spectral variable t s ,z in CPM and six-vertex model, respectively, we denote 



i-t 



H(t) 



t ■= r]- l t s , in CPM, 



i-t_ 

1 — t, t := qz, in six — vertex model. 



For non-zero parameters v±, ■ ■ ■ , v m , distinct and not containing a complete iV-cyclic string (which 
means {lu^v} j e 2, N ^ or some v / 0), we consider the non-zero polynomial 

m 

F(t)=n(i+u J —v i t), 
i=i 

where no is the integer defined by no := —2, 1 for CPM and the six- vertex model, respectively. 
Note that as before, u = q 2 in the six- vertex model case. For an integer r with < r < N — 1, and 
the polynomial F(t), we define the function 

P(t) is invariant when changing t to uit, hence depends only on t N . The polynomial criterion for 
P(t) is that the v^s in F(t) satisfy the relations 

—, 20 \' „) m+i -i ■ * = !,•••, m. (3.6) 

H{-uj 2 Vi ) L F(—u> 2 +L v i ) 



2 The Q-operator here differs from the Q cp in [28] by a scale factor: (77 2 :r s ) p °(?7 2 y a ) p <> /_i 3 >J Q(s) = Q cp (s). 
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Indeed the above equation is the Bethe equation for the superintegrable CPM ((4.4) in [1], (6.22) 
in [7]) and for the six- vertex model where r = ^ — m (= \S Z \) (mod N) with < r < N — 1. 
Then the polynomial P(t) in (3.5) is the evaluation polynomial for the symmetry algebra of the 
model [18, 28] 3 . 

Based on the common feature of the Bethe equation (3.6), and known functional relations in 
CPM about operators T^\t), Q(s) in (3.4) (3.3) and symmetries U, C in (3.2), we now derive func- 
tional relations for both CPM and the root-of-unity six- vertex model. The (t) and symmetries 
U, C in six- vertex model are defined by 

TW(t)= ziT(z), (t:=qz), ^ 
U : z i— > uz, C ' s ' ^ —s, (s := zs), 

where T(z) is the transfer matrix (2.2). Note that T( 2 )(i) is a ^-polynomial operator, and in terms 
of the variable s, U is defined by U(s) = qs. 

By the Bethe equation (3.6), the T( 2 )-eigenpolynomial, denoted again by T^ 2 \t) (when no 
confusion could arise), satisfies the relation 

T^{t)F{t) = uj- r H(t) L F{uj^t) + H(u- 1: t) L F(ujt). (3.8) 

Define the T^-operators for non-negative j recursively by setting T^°\t) = 0, T^\t) = i?(oj _1 t) i , 
and the fusion relation 

r^T^V'- 1 ?) = ^HiJ'H^T^it) + H(J- 2 t) L T (j+1 \t), j > 1. (3.9) 

Express T^ +1 ^ in and T^ -1 ) from the above relations; then by induction argument, one obtains 
the eigenpolynomial expression of T^\t): 

T«\t) = F{u-H)F{ui-H) £ g U k Q " , j>1. (3.10) 

By this, follows the W-periodic' relation of T^'s: 

v ) F(t) v ; w F(t) 

in particular, the case j = 1 yields the boundary fusion relation 

= uj~ r T^- l \u:t) + 2H{u- l t) L . (3.11) 

Note that the fusion relations, (3.9) and (3.11), are equivalent to those for t^'s in the superinte- 
grable CPM case ((53) in [28]). 

3 The evaluation function of six-vertex model in [18] is given by (3.9) with Y(v) = ^^Sq^v + 2(j + l)i7), 
where a(v) = ™ — 2 — ' 7 .^ in (2.47). In terms of variables z.t and Bethe roots ViS here, a(v) = 

2 2m - L (nr =1 ^)(g- i ^^ +m TT ^ -rj^gffe^) ' which implies Y(v) = 22m_i (n™ i ^ +m+r m 
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We now discuss the Q-operator, i.e. a family of non-degenerated commuting operators Q(s) 
depending on the variable s (which is algebraically related to t by (3.1) or (3.7)), such that 
[T(t), Q(s)] = with the following TQ-relation: 

(t)Q(s) = uj^'Hit^QiU^s) + H{u~ l t) L Q{Us), (3.12) 

where U is the s-symmetry in (3.2) or (3.7). Note that by comparing (3.8) with (3.12), the t- 
polynomial F(t) is expected to be a factor of an eigenvalue of Q{s). Using (3.12) and (3.9), one 
finds the following T^Q-relation by induction argument for j > with = 0, 

T ( i){t) = Q(U- 1 s)Q(U j - 1 s) 3 J2 (^ kr H(io k - 1 t) L Q(U k - 1 s)- 1 Q(U k s)-A. (3.13) 

fc=0 ^ ' 

Indeed for a Q-operator in (3.12), the T^> 'Q-relation (3.13) is equivalent to the fusion relations, 
(3.9) and (3.11). The TQ- and T^Q-relations here are equivalent to the r^ 2 ^Q cp - and T^Q cp - 
relations in CPM case ((54) (55) in [28]). In the case of six- vertex model, the TQ-relation (3.12) is 
valid for any Q-operator (up to a suitable normalizing factor), so is the T^Q-relation (3.13). 

3.2 The QQ-relation and Q-functional relation 

Now we discuss the constraint on the Q-operator to encode the essential features about the sym- 
metry of CPM and the root-of-unity six- vertex model. First we note that by (3.10), the T^)- 
eigenpolynomials and evaluation polynomial P(t) in (3.5) are related by 

T^(t) + u>- jr T( N -i\u j t) = (w _1 t) r J P(o;- 1 t)F(^- 1 i)P(a;- 1 F), < j < N, 

T( N \t) = {ui-HyFiu-HfPiu-H), (3-1 ' 

which in CPM case are equivalent to the relations between r^-ei genvalues and the evaluation 
polynomial ((73) in [28]). We expect the Q-operator in accord with (3.14) and the expression of 
P{t) in (3.5), which leads to the following constraint of Q(s): 

,JV-1 H(m k - 1 T) L (m k t)- 



— r 



T (N) /I\ _ 7r n( jj-l q) 2 WV-1 H(^-H) L {uH)-r 
1^ >[t)-ttJ{U S) ^fc =0 Q{U^s)Q{U k s 



The [/-invariant operator J2^=o H Q(u k -%)Q(U k s) ^ as ^ ne eigenvalue uj~ r P(uj~ 1 t). We now use the 
symmetry C in (3.2) or (3.7) to state the functional relation required for Q(s): 

N-l 

P(t )F(t ) 2 Q(Cs) = Q(s) Y, H(co k ~ 1 t) L (u; k t)- r Q(U k - 1 s)- 1 Q(U k s)- 1 , (3.15) 

fc=0 

with to = 1)0 in CPM and the six-vertex model, respectively. Here, we abuse the notation by 
using P(to)F(to) 2 to denote the diagonal operator acting as the (P(to)F(to) 2 )-m\iltiple on the 
-eigenspace with the -eigenvalue determined by F(t). Then, the above constraint relations 
of and Q yield the following QQ-relation: 

Pit^Fitoftyiu-^Mw^Cs) = T0)(t) + w-^t^-j V'*)> o < j < n, 

P(t )F(t ) 2 t r Q(U~ 1 s)Q(U~ 1 C S )=T( N \t). [ ' ' 
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In the case of CPM, (3.15) and (3.16) are valid, and equivalent to the Q cp -functional and Q cp Q C p- 
relation in [29] (for formulas (33) (56) by using the equality N L u Pb e iP = P(1)F(1) 2 ). Indeed, 
relations in (3.15) and (3.16) are all equivalent by the following result. 

Theorem 3.1 For an integer j between and N, the jth QQ-relation in (3.16) is equivalent to 
the Q-functional relation (3.15), hence equivalent to the whole set of relations (3.16). 

Proof. By (3.13), one has 

V7) = o^-MQdrMES- 1 

By substituting T^>) of (3.13) and the above expression into (3.16), follows the equivalence of the 
jth relation in (3.16) and the Q-functional relation (3.15). □ 



For the rest of this paper, we consider only the root-of-unity six-vertex model in sectors S z = 
(mod N), where = r = ^ — m. The polynomial F(t) and Bethe equation (3.6) are now expressed 
by 

m / i —\\L m — 1 

F(t) = U(l + v t z), (11 =-Y[ Vl UJ - for z = l,...,m, (3.17) 

with t = qz. Hence, F(0) = 1, and P(0) = N. Note that the equation (3.17) is unaltered when 
replacing v { by v~ l . Denote F^(t) = Ui^ii 1 + v^z). Then the polynomials P(t),P^(t) of (3.5) for 
F(t), F^(t) resp. satisfy the following 'reciprocal' relation: 



t L - 2m P{t- 1 ) = ([[v l )- 2 P\t). 



i=l 



By this and P^{0) / 0, the t-degree of the polynomial P(t) is given by 

degP(t) = L — 2m. (3.18) 

The -eigenvalue is determined by the Bethe polynomial F(t) via (3.8), now in the form 

T^(t)F(t) = (1 - t) L F{oj- l t) + (1 - u-H^F^t). (3.19) 

The non-degeneracy ofT( 2 )(t), i.e. a T^ 2 ) -eigenvalue with one-dimensional eigenspace, is equivalent 
to the constant function for the polynomial P{t) associated to F{t). As a consequence of (3.18), 
we have the following characterization of non-degenerated -eigenvalues: 

Lemma 3.1 Let T^ 2 \t) be a -eigenvalue corresponding to a Bethe polynomial F(t) in (3.17). 
Then T( 2 )(t) is non- degenerated if and only if the degree F(t) is equal to 

□ 

Now the r( 2 )Q-, QQ- and Q-functional relations, (3.12) (3.16) and (3.15), take the form: 

TW(t)Q(s) = (1 - t) L Q(q- 1 s) + (1 - uo- l t) L Q{qs), s := z\, (3.20) 
NQ{q- 1 s)Q(-q^~ 1 s) = T^\t) + T^ N ~^\uH), < j < N, (3.21) 

NQ(-s) = QOOEj^l - u k ^t) L Q{q k ^s)^Q{q k s)-\ (3.22) 

In the next section, we produce a Q-operator with the above properties. 
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4 The Q-operator for Six- vertex Model at Roots of Unity 



We first construct the Q-operator of the root-of- unity six- vertex model in subsection 4.1 by following 
Baxter's method of producing the Q-operator of the eight-vertex model in Appendix C of [3]. 
Then in subsection 4.2 we verify the Q-functional equation (3.22) for the size L < 6 by direct 
computations. Finally, in subsection 4.3 we show the functional relations hold for a general L by 
using the explicit forms of the Q-operator and fusion matrices. 

4.1 Construction of the Q-operator 

The L-operator (2.1) is the matrix with the C 2 -auxiliary space and the following (^-(quantum- 
space) operator-entries: 



where 



06/ ' \ c ' \ / ' \ a 



i i -i -i i i ii -i -i i 

a(z2) = z?q 2 — z z qz , b = b(z2 ) = z'2 q2 — z 2 q 2 ; c = q — q . (4.1) 



Since the chain-size L is even, the entries of the transfer matrix T(z) are z-functions. Hence T(z) is 
unaltered by changing the variable to — z? in the L-operator, i.e., T(z) = ti a ux ®f=i z?). 4 
Consider an S-operator with C^-auxiliary space and C 2 -quantum space, S = {Sij) i j & z N > sucn 
that the C 2 -operators Sjj are of the form 

SiJ = VijTij, ViJ = ^ ^ ^ ' T M = ( T 1' T 2)> ( 4 - 2 ) 

where Vj's depend on the z- variable, and Tj's are constants in C. The general form of the matrix 
Qr we shall use is 

L 

Qr = tr qn S^) , = S at the site £, 
e=i 

with TQr = tr q2 ^ qn (®f =1 \Je) and = U at the site £, where U is the following matrix 
associated to the L- and S-operator with C 2 <g> C^-auxiliary space and C 2 -quantum space: 

U _ I ^o,oS L 0i iS \ 
\ LifiS Li ; iS J 

The operator TQr will decompose into the sum of two matrices if we can find a 2N by 2N scalar 
matrix M (independent of z) such that 

M~ 1 UM = ( A B V 
V D } 



4 The functions, a(-z^), b(-z^) and c, correspond to Boltzmann weights of the six vertex model in the face-model 
description in [9] (22): W 6v {a, /3, 7, S) = 8{a,- 1 )(zq) 1 a(-z?) - 8{f3, S)(zq- 1 ) a b(-z^). 
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As in [3], the above required form is unaffected by postmultiplying M by a upper blocktriangular 
matrix; together with a similar transformation of S, we can in general choose 



M 



In 
5 I 



, <5 = dia[<5o,-",<Sjv-i]- 



N 



Hence, 

M" 1 UM=f Lo,oS + L ,iS<5, L 0)1 S \ 

\ -5L 0t0 S + Li ;0 S - 5L ^SS + Li,iS£, -<5L ,iS + Li,iS J 

and the condition on Sjj's in (4.2) for the vanishing lower blocktriangular matrix is 



-aSi + b8j c 
—c5i5j a5j — b5i 



Vl 
V2 



(ri,r 2 ) = 0. 



For non-zero Sjj, this in turn yields 



= (a 2 + b 2 - c 2 )SiSj - ab(5 2 + 5 2 ) = (z + z^ 1 - (q + g~ 1 ))(5 l - q5 j )(S l - q^S 



'3)1 



hence Si = q ±x Sj and 



c : (aq — b)5j if Si = qSj, 
c : (aq^ 1 — b)Sj if Si = q~ 1 Sj- 



We choose the L-operator to be L{— za ) in the above argument without affecting T(z), i.e., replacing 
a, b by a(— z*), b(—Z2), and then set 



5 = dia[l, q, ■ ■ ■ , q N ~ l ], i.e. Sj = q 3 for j = 0, . . . TV - 1. 



(4.3) 



Then Sij is equal to zero except j — i = ±1 G Zat, in which case the two- vector ratio for Vjj in 

j — i i— j 

(4.2) is given by v\ : t>2 = 1 : z 2 g 2 Sj. Hence, we can set 



((zq) 452, ( Z q) ^ q 2 y if j - i = ±1, 



with the S-operator in the form 

/ 

Si,o 

s 







So,i 

Si )2 











otherwise, 



(4.4) 



-1 \ 





SjV-2,JV-l 





V SjV-1,0 • • • SjV-l,iV-2 

One can show the (i,j)-th entry of the diagonal block matrices of M _1 UM satisfies the relations 



1 \ i- 



(L 0fi S)ij(z) + (L 0jl S)ij(z)Sj = a(-Z2)q 2 S i:j (zuj), 
(-5L 0A S)ij(z) + (L 1A S)ij(z) = bi-z^q^T^Sijizu;- 1 ), 
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which, by L even, imply 

T(z)Q R (z) = b(z^) L Q R (zu- 1 ) + a(z^) L Q R (zu;). (4.5) 

We now replace Qr by Ql, and form the product QlT by using Ql = tr^~,iv(0^ =1 S^) where the 
operator S = j e % N with Sjj in the form 

Si j = njy/ij, %j = ( f 2 ) ' = ^ %U%2) ' ^ 4 ' 6 ^ 

We repeat the above working, replacing S, L(-z^) by S, L(z?), but using the same M with 5 in 
(4.3). We find that 

~ . = | ({zqY-r-qi^zq^q^) if j - i = ±1, 
\ otherwise, 

and 

Q L (z)T(z) = 6(z5) i Q L ( zw - 1 ) + a(^) L Q L (W). (4.8) 

Note that the z-dependence of Qr(z),Ql(z) in (4.5), (4.8) indeed means operator-functions of z*. 
In order to construct a Q-operator commuting with T, as in [3] (C28) it suffices to show 

Ql(w)Qr(z) = Q l (z)Qr(w). (4.9) 

To do this we follows the method in [9] (formulas (48) (49)) by considering the product functions 
of (4.4) and (4.7), f(w,z\i,j;k,l) = \j i:j (w)\i ki i(z): 



f(w,z\i,j;k,l) 



(wq) 3 ± (zq) 4 q 3 2 + (wq) ± (zq) t q ^ if \j — i\ = \l — k\ = 1, 



otherwise, 

and looking for an auxiliary function P(w, z\n) for n £ Z^r such that 

f(z,w\i,j;k,l) = P(w,z\k-i)f(w,z\i,j;k,l)P(w,z\l - j)' 1 . (4.10) 

Since the product Ql(z)Qr(w) differs only the boundary contribution when interchanging w by v', 

Ql(z)Qr(w) = P(w, z\ki - ii)Ql{w)Qr(z)P(w, z\k L +i - zl+i) -1 , 

this in turn implies the commutation relation (4.9) as the P-factors cancel out by the periodicity 
of boundary condition. There are four cases, \j — i\ = \l — k\ = 1, to consider for the above function 
P. The relation (4.10) automatically holds for j — i = I — k, and the other two cases yield just one 
condition for P: 

P(w, z\n + 1) (wq) 2 + q n (zq) 5 
P(w,z\n-1) ( Z q)^ +q n (wq)^' 
By ./V odd and q being an Nth root of unity, the function P is determined by the relation 

P(w,z\2n) -fj (wq)^ +q k (zqfi 

—fT, uyT = I I 1 ; r tor n G L N- 

P(w,z\0) ^(zq)2 +qk(wq)2 
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By this, (4.9) holds. If we define 

Q(z) := Q R (z)Q R {z )- 1 = Ql(z )- 1 Ql(z), (4.11) 

then 

[T(z),T(w)] = [T(z),QH] = [Q(,),QH=0, (4.12) 

and 

T(z)Q(z) = biz^Qizcj- 1 ) + a{z^) L Q(zuj). (4.13) 

— c z — <j z 

Note that Vij, Vjj in (4.4), (4.7) satisfy the relation, q~^~\ii j = Vj+ij+i, Vjj = Vj+ij+ig - ^ - . It is 

a z —<j z 

convenient to choose Tjj, %j in (4.2), (4.6) with the relations, q~Tij = Tj+ij+i, q^~%j = fj+ij+i. 
Then, 

-° z c si c £±? ? 

Q 2 2 - =>i+l,j+l, Q 2 3i,jQ 2 - ^i+lj+l- 

So Tjj's are determined by to,i, t^jv-i; the same for Tjj's by ?o ; i, t^jv-i- Since (4.11) is unaffected 
by post- and pre-multiplying Qr(z), Ql(z) respectively by constant matrices, we may set 

T~0,1 = (1|, To,JV-l = ( — 1|, %,1 = To,JV-l = | ~ 1), 

where (±1| are the dual bases of | ± 1). Hence nj's and %,/s are zeros except j — i = ±l (mod N), 
in which cases one has the expression 



i(j-i) , , ^ -'(j-Q 



'2J 



2 



^j=9 2 b" — iorj-i = ±\. 



Then by (4.4) and (4.7), the only non-zero Sjj , Sjj are 

Sij{z) =q^E a =±iq^ 1 (zq)^~ jl \a}(j iorj-i = ±l (mod AT), 
Sij(«) = 9^~ li E / 3=±i9^(^) _LLrlM |i -*)</?|, forj-i = ±l (mod AT). 

Note that up to (— l)-power factors, the entries of Sjj(z), Sjj(z) are the same as the weights in 
formulas (27) (43) of [9]: indeed by z = e v ,q = e x , and the identification a,j,i with b — a,c,d, 
Sij(z) corresponds to Wq(v\a,b,c,d) in [9]; the same for Sij(z) and Wq(v\a,b,c,d) by identifying 
i,j, (3 with a,b,c — d. 

We now calculate the matrix form of Q/j and Q^. The relations, (4.14) and Qr — ®^ = iSi^^ +1 ( 
with ^ e Ztv and = i\, in turn yield 

Qr = 2^ i 2 i 2 4 ^ |"^) ^ (v+i - vl, 

if a t =±l 

where the prime-summation means the sum of those indices it's with — ii = ±1. Denote 
Pi = ii+i — i( = ±1 for 1 < i < L, and % = i\ G Zjv- Then E^ = i^+i — *i = (mod A 7 "), and 
Qr, as an operator of V in (2.3), is expressed by 

Qi? =Ea £ ,ft=±i9 2 4 Ef=o 9 2 ®/ l<*<} ®< </%|, 

= iv Ea^A=±i9 2 ( z <?) 4 ®e\(*i) ®e(Pe\- 
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Since Ei< s <kl PsPt = J2s Ps J2e >s Pt = -Y. a Ps{Yz< s Pi) = -L- J2i<e< s <L PtPs (mod N), one 
has 

<-! -T 

EAE^ E T (mod TV), (4.15) 

£ s=l 

by which we obtain the matrix expression of Qr: 

Qr = Nq— 1 ^ N) 3 ^ K> ®* (/H ( 4 - 16 ) 

where indices o^,/% for 1 < I < L in the summation are ±1 with J2t a e = HtPt = (mod AT). 
By L even, the above expression of Q# indeed defines an operator with s-polynomial entries, where 
as before s = z?; and we shall also write Qr = Qr(s). We repeat the above calculation, replacing 
Qr, Sjj by Ql, Sjj-, then we find the matrix form of the F-operator Q^: 

l ~ 2~^£ a ^ Yl —\ @ s ~ a e@£ 

Q L (s) = Nq~ V 5 ~ J= ( z l) 3 ®z \ a i) ®t {Pel ( 4 - 17 ) 

Hence Ql and Qr differ by a scalar, both commuting with the spin-inversion operator R, 

Ql(s) = q%Q R {s), [R, Q R (s)} = [R, Q L (s)} = 0. (4.18) 
Using J2e a i = J2e Pi = ® (mod AT), we find that the spatial translation commutes with Qr, Ql: 

[Sr,Qr(s)] = [S r ,Ql(s)]=0. (4.19) 

4.2 Justification of the Q-functional relation for L = 2,4, 6 by direct computa- 
tions 

Here, we demonstrate by examples that the Q-operator in (4.11) satisfies the QQ- and Q-functional 
relation in subsection 3.2. We modify Q^, Ql by a normalized factor in accordance to the six-vertex 
transfer matrix in (3.7), and write 

Qr{s) = z^Qr(s), Q l {s) = z^Q l {s). 

By (4.18) and (4.19), the s-polynomial operators Qr,Ql commuting with R,Sr with Ql{s) = 
Q^Qr( s )- Denote 

Q(s) := Qr(s)Q r (0)- 1 = Ql^QUs). 

As in (4.12) and (4.13), Q(s) is the Q-operator associated to so that the T^Q-relation (3.20) 
holds. By (4.16), Qr(0) = Nq~ld, and the matrix expression of Q(s) is given by 

Q(s)= 2^ q 2 q 4 s 2 ®e\a e )®e{Pt\- ( 4 -20) 

a e ,f3 e =±l 

By (4.15), Q(s) is a s-polynomial operator of degree L with Q(0) = Id and the leading coefficient = 
the spin-inversion operator R. Hence [R, Q(s)] = 0. Furthermore one can show Q(—s) = SQ(s)S 
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where S = <S>£ cf ■ Indeed with Tjj in (4.2), we define Tjj = ±Tjj according to j — i = ±1 (mod N). 
Replacing by Tjj in the definition of Qr in subsection 4.1, we obtain the operator Qr, hence 
Qr(s) ■= s~Qr(s). The operators Qr and Qr are related by Qr(—s) = SQr(s), which implies 

Q(s) = Qd-s^Rioy 1 = SQr(s)Qr(0)^ 1 Qr(0)Qr(0)^ 1 = SQ(s)Q R (0)QR(oy 1 . 

Setting s = in the above equality, one obtains Qi?(0)Qi?(0) -1 = S, therefore Q(— s) = SQ(s)S. 
Note that as an operator of V, one can write the Q(s) in (4.20) in the form 



Q(s) 



L 
S2 



E 

at,/3e=±l 



<3e \ae) ®£ {Pe\- 



Up to s~-factor, the above operator coincides with the six-vertex Q-operator (101) in [4] which 
holds for a generic q with J2e a z = J2ePe = 0- Hence, (4.20) can be considered as the Baxter's 
classical six-vertex Q-operator in the root-of-unity case. 

We are going to show the above Q-operator satisfies one, hence all by Theorem 3.1, of relations 
in (3.21) and (3.22) for cases of L < 6 by computational methods. For simple notation, we shall 
also write basis elements of the vector space V in (2.3) by 



(4.21) 



First, we consider the case L = 2, where V is the 2-dimensional space with the basis, \a, —a) 
(a = ±1). Then T(z) in (2.2), and Q are expressed by 



T(z) = 

r( 2 )(i) 



{z + z~ l ){q b +q~ b )-2(q 2 <g> q 2 + q 2 <g> q 2 ) + ( q - q- l ) 2 (a_ 
2(z 2 + l-z(q + r 1 ))(|l, -1| + I - 1, !)("!, 1| 

l-afi 



cr+ + 0+ <g> (j_), 



+z{q - q- 1 ) 2 ^ -1,1)(1,-1|+|1, 1 1|), 



Q(s) = E a , =±iz—\a,-a)((3,-p\. 

Note that Q(s) is a z-polynomial operator. The eigenvectors of Q and are: |1, —1) + | — 1, 1) 
, |1, —1) — I — 1, 1), with the Q- eigenvalue (1 + z) ,(1 — z), and T( 2 )-eigenvalue 2(z 2 + 1 — z(q + 
q -1 )) + z(g — q^ 1 ) 2 , 2(z 2 + 1 — z(g + q" 1 )) — z(q — q" 1 ) 2 , respectively. The above Q-eigenvalues are 
indeed the polynomial F(t) for the Bethe equation (3.17), in which case the Q-functional equation 
(3.22), equivalent to (3.5), automatically holds. 

For L = 4, V is of dimension 6 with the basis, vj^, 1 < j < k < 4, where vj^ = ®\ott) with 
ct£ = 1 for £ = j, k, and —1 otherwise. Denote Vj~ k = Vj : k + Rvj : k, v~ k = Vj^ — Rvj t k- The matrix 
form of Q for the basis {^2, v 3A ,v 1:3 ,v 2A , v lj4 , v 2 , 3 } is 

/ 1 z 2 

z 2 1 

z z 1 z" 

z z z 2 1 



Q(s) 



z z q 
J2 



q l z 



qz 



z z 



qz 


q l z 






~ l z 


qz 


z 


z 


z 


z 


1 


z 2 


z 2 


1 
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Hence with respect to the basis {v^ 2 j , 3 j , 4 j , 2 j , 3 j , 4 j }, Q is expressed by 



Q00 = 



V 



l + z 2 
2z 







2z 
1 + z 2 
2z 






(q + q Y )z 
2z 
l + z 2 






(3 







l-z 2 


1 -q)z 








l-z 2 




o 




<r 
o 



\ 



and the transfer matrix T(z) is given by 



T(z) 



V 



2a 2 b 2 
2abc 2 
(a 2 c 2 + b 2 c 2 ) 









2abc 2 
{2a 2 b 2 + c 4 ) 
2abc 2 









( a 2 c 2 + 6 2 c 2) 

2a6c 2 
2a 2 6 2 













2a 2 6 2 


(b 2 c 2 



(2a 2 b 



l-z 2 








2 1.2 



(4.22) 









(a 2 c 2 



a 2 c 2 ) 



b 2 c 2 ) 


2a 2 b 2 



) 



where a,b,c are in (4.1). The eigenvectors for and Q are: 



v l,2 



v iA> v t,2 + v iA ~ z(« + 



g 1 ± \/ q 2 + q 2 + 34)^^, ^ 2 ± iu 14 ,v 13 , with distinct eigenvalues. Indeed the Q-eigenvalues, 
1 + z 2 - (q + (T 1 ^ , 1 + z 2 + 9±^rVp^±34 Z) 1 - z 2 ± i(g - q -^) z ,1-z 2 , are F-polynomials 
for the Bethe equation (3.17) with m = 2, (true as well for an arbitrary q which is not necessary a 
root of unity). Then the z-polynomial expression of Q yields the Q-functional equation (3.22). In 
fact, one can directly verify (3.22) by using the matrix form of Q(z). For example for N = 3, the 
Q-functional equation (3.22) of size L is equivalent to the following one: 



3Q(-s)Q(qs)Q(q 2 s) = (1 - z) L Q(s) + (1 - qz) L Q(q 2 s) + (1 - q 2 z) L Q(qs). 



(4.23) 



Then the matrix (4.22) satisfies the above relation for L = 4. 

In the cases L = 2,4, all -eigenvalues are non-degenerated, and the Q-eigenvalues are the 
z-polynomial Bethe solutions of (3.17). This kind of relation between Q-eigenvalues and Bethe 
solutions indeed holds for all non-degenerated -eigenvalues for a general L by assuming Bethe 
ansatz of the six-vertex model. 

Theorem 4.1 LetT^ 2 \t) be a non- degenerated -eigenvalue for an arbitrary (even) L. Then 
its eigenspace is spanned by a vector in V with S z = 0, and the corresponding Q-eigenvalue is a 
z-polynomial which defines a Bethe solution F(t) of (3.17) with m = ^. As a consequence, the 
Q-eigenvalue satisfies the Q-functional relation (3.22). 

Proof. By Bethe ansatz, there exists a Bethe polynomial F(t) for a solution of (3.17) with 

m = ^. This im- 



m 



< h so that T^ 2 \t) and F(t) satisfy the relation (3.19). By Lemma 3.1, 



plies the one-dimensional eigenspace with the eigenvalue T^ 2 \z) has the basis element v with 
\S Z \ = ^ — m = 0. Then Q(s)v = Xv. Express v as a linear combination of the standard basis, 
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v = J2/3 1 +...+i3 L =o a i3u-,i3 L \Pi, ■ ■ ■ We may assume one of a /3u ...,p L , s taking the value 1, say 

°7i,---,7l = 1' hence by Rv = ±v, a_ 7lv .. 5 _ 7i = ±1. The relation (4.20) yields 

A = \7i) ■ ■ ■ > Jl\Q{ s )v = }^ a Pu-,PL<l 2 Q 4 z 4 

/3i+...+/3z=0 

In the above sum, ^i^ 7 ^ are integers between and ^, taking the value ( \) only when 
0t = -Yt( ~lt resp.) for all £. By (4.15), E^EfciT* = (mod TV), hence A (= \{z)) is 

a z-polynomial of degree ^ with A(0) = 1 and the leading coefficient = ±1. Then the equation 

(3.19) holds with F(t) := X(z). Claim. The roots of the polynomial X(z) do not contain a complete 
iV-cyclic string, consequently, follows the result. Otherwise, by factoring out those roots of X(z) 
which form complete TV-cyclic strings, one obtains a polynomial F'(t) with degree less than ^, and 
roots containing no complete TV-cyclic string, such that T 2 \t) and F'(t) again satisfy (3.19). Then 
F'(t) is a Bethe polynomial for a solution of (3.17) with m < ^, which by Lemma 3.1 contradicts 
the non-degeneracy of TW(t). a 

As a corollary of the above theorem, the Q-functional relation holds when all T( 2 )-eigenvalues 
are non-degenerated. 

Proposition 4.1 The necessary and sufficient condition for the T^ -matrix having no degenerated 
eigenvalue is the relation ^ < TV. In this situation, the Q-operator (4.20) satisfies the Q-functional 
relation (3.22). 

Proof. The condition ^ < TV is equivalent to V being a subspace with the spin S z = 0, hence by 
Lemma 3.1, equivalent to all -eigenvalues being non-degenerated. Then the results follow from 
Theorem 4.1. □ 

Remark. The conclusion in the above proposition about non-degenerated T-eigenvalues and the 

corresponding eigenvalues of Q holds for all Q-matrices, not a special property of the Q-operator 

(4.20) . The computational argument we present here has further enhanced the correct form about 
the operator (4.20). □ 

We now discuss the case L = 6. For TV > 5, the Q-functional relation holds by Proposition 4.1. 
So we need only to consider the case TV = 3, where V of dimension 22 with the basis consisting of 
v := <8>|l),v' := <8>| — 1) Vi,j,k, l<i<j<k<6 (:= the vector in (4.21) with only the k-th 
factors being |1)), and the dual base, v*,v'*, v* • fe s. 

By Theorem 4.1, there is only one degenerated r^ 2 ^-eigenvalue, given by the 'pseudo-vacuum' 
v with the eigenvalue (1 — uj^t)^ + (1 — i) 6 . Denote the corresponding eigenspace by W, which is 
invariant under the Q-operator. So it remains to show that functional equation (4.23) holds for Q 
when restricting on W. With F(t) = 1 in (3.5), one has the evaluation polynomial 

P(f) = l-20P + F\ (4.24) 
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with the Drinfeld polynomial, 1 — 20x + x 2 (see also (8.25) in [15]). Denote the 5^-invariant vectors 
in V, 

5 5 5 

UO := Ui )3 ,5 + ^2,4,6, «1 := V lk,2 k ,3 k , ™ •= Yl U lfc,2fc,4 fcJ «>' == £ «l fc ,2 fc ,5 fc , 

fc=0 fc=0 fc=0 

where = 6 k (i) with # the cyclic permutation of {1, • • • ,6} sending k to k + 1 by the identification 
7 = 1. Then i?(uj) = Ui (i = 0, 1), !?(«;) = to'. The s/2-loop-algebra structure of the space W is 
determined by the values of the generators (2.4) on v,v', expressed by 

S-W(v)(= T +( - 3 \v')) =u + Ul +qw + q 2 w', T'^(v){= S +(3) (v')) = u + «i +q 2 w + qw', (4.25) 

with the relations (S-^) 2 {v) = 20v', (T-^) 2 {v) = 20v', T-^S'^iv) = S'^T'^iv) = 2v'. 
So W is the 4-dimensional vector space with {v,T~( 3 \v), (v), v'} as a basis. We are going to 
determine the matrix form of Q with respect to this basis. By (4.20), for any basis element b of V 
with b* as its dual element, 

<6*|Q(s)|6) = 1, (b*\Q(s)\(Rb)) = z 3 , 
and the following relations hold: 

(ai,..., a 6 \Q(s)\v) = q 2 z * =z 2 (v'\Q(s)\ai, . . . , a 6 ), 



Hence, 



{a!,. . . ,a 6 \Q(s)\v') =q 2 z a = z 2 { v \Q(s)\ai, . . . , a 6 ). 



(v*\Q(s)\v) = {v'*\Q(s)\v') = 1, W*\Q(8)\v) = {v*\Q(s)\v') = z\ 

(v*i, 3 , 5 \Q(s)\v) = (v* lt3t5 \Q(s)\v>) = z§, (vt X3 \Q(s)\v) = (v* 1X3 \Q(s)\v>) = z§, 
(vl 2A \Q(s)\v) = K i2)5 |Q(a)K> = <M « 2)5 |Q(s)b> = « 2 , 4 |Q(s)K> = • 



Since [Q,Sr] = 0, and v,v' are 5_R-invariant vectors, the above relations yield 



Q(s)v = v + z 3 v' + z2uq + 22^2 + q 2 z?w + qz^w' = v + z?T ( 3 \v) + z 3 v', 
Q(s)v = z 3 v + v' + Z2ti + z2^i + gz2tt; + q 2 z2w' = z 3 t> + 212 £ - ( 3 )(i;) + v' . 



(4.26) 



By 

{^!,2,3l<9( S )l^l,3,5> = «2,4l<2(s)K3,5> = {v*2, 3fi \Q(s)\ Vl ^ b ) = Z, 
«2,3l<5( s )l w 2,4,6) = K,2,5lQ( s ) 1^1,3,5) = (^2,3,5 1 <2( s ) l«l,3,5> = Z 2 , 

the S'ij-symmetry property of Q yields 

Q(s){u ) = 2z*v + 2ziv' + (1 + z 3 )u + {z + z 2 )( Ul +w + w'). (4.27) 

By 

«3,5l<2(s)K2,3) = (V1,2A\Q{S)\V1,2,3} = Z, (^,4,6 \Q(s) |«i,2, 3 ) = «4,5 1 <9( S ) h,2,3> = Z 2 , 
{vl 3A \Q(s)\v 1;2 ,3) = <«l,3,6lQ(*)l«l,2,3> = 9*, (^,2,6l<9( S )kl,2,3) = (^2,3,5 1 <9( S ) |wi,2,3> = ? 2 «, 
(^3,4,5lQ( S )l U l,2,3) = <«2,5,6lQ( s )l"l,2,3) = 9* 2 , «5, 6 IQ( s )kl,2,3) = K 4 , 6 IQ( S )K2,3> = 
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the R, ^-symmetries of Q yield 



Q(s)(ui) = Qz^v + 6z%v' + 3(z + z 2 )u + (1 - z - z 2 + z 3 )m. 



(4.28) 



Similarly, using 

«3,5lQ( s )K2,4> 



«5,6lQ( s )l v l,2,4) = (W3,4,6IQ( S )I U 1,2,4) = (l>2,5,6 1 Q( S ) I V l,2,4> = ( V 2,3,G 1 I w l,2,4> 

K 4)6 |Q(s)K2,4> = K2,3lQ( S )K2,4> = (u£ 4)5 1 Q(s) | t>l,2,4> = K.2,5 1 <2( S ) K2,4> = K j3)4 | Q(s) \ V 1>2 ,4,) 



(v*2,4,5\Q(s)\v1,2a) = Q 2 Z, <«l,2,6lQ(*)l«l,2,4) = K 4 ,6 1 <2( S ) |«1,2, 



gz, 



K3,5lQ( S )l^l,2,4) = Q 2 Z 2 , {vt t5fi \Q(s)\v lt2 ,A) = «3,6 1 <9(«) |wi,2,4> = 9^, 



(4.29) 



(^3 )4 |Q(s)|wi i2 ,4) 
K4,5lQ( S )l U l,2,4) 

one obtains 

Q(s)(w) = 6qz^v + 6q 2 z%v' + 3(z + z 2 )u + (1 + z + z 2 )w + (z + z 2 + z 3 )w' , 
Q(s)(w') = 6q 2 z^v + 6qz%v' + 3(z + z 2 )u + (z + z 2 + z 3 )u> + (1 + z + z 2 )w' . 

By (4.25), (4.27), (4.28), and (4.29), we have 

Q(s)(T~ 3 (v)) = 20z%v + 2ziv' + (1 + z 3 )(u + ui) + (q 2 + <?z 3 )u? + (q + g 2 z 3 )u>', 

= 20ztw + 2z§ u' + T~ 3 (v) + z 3 ^- 3 ^), 
Q(s)(S -3 (t;)) = 2zlw + 20z^v' + (1 + z 3 )(u + ui) + (q + q 2 z 3 )w + (q 2 + qz 3 )w' 

= 2ziv + 20zf v' + z 3 T~ 3 {v) + S _3 (v). 

Hence by (4.26) and (4.30), with respect to the basis {v, T~( 3 ) (v), S~( 3 ^ (v), v'} of W the matrix 
form of Q is 



(4.30) 



Q(s) 



3 

Z2 





20z§ 2z§ 



1 



,3 \ 





3 
Z2 



V 



20zi 



i 

Z2 . 



Z2 

Z5 

3 
2 

3 



z" 2z2 

with the following eigenvectors and eigenvalues: 

Eigenvector, Eigenvalue 

3V2(v-v')+T-W(v)- S-W(v), 1 + 3V2. " 

3V2(v - v') - T-^{v) + S'^iv), 1 - 3^2, 

V22(v + v') +T~( 3 \v) + S~ (3 \v), 1 + V22zl + 

^{22)(v + v') - T~^(v) - S-^iv), 1 - v^zf + z 3 . 

Then the above Q-matrix satisfies the Q-functional equation (4.23) for L = 6. The evaluation 
polynomial (4.24) is related to the Q-eigenvalues by the factorizations: 

P(t) = (1 + 3y/2z* - z 3 )(l - 3\/2zt - z 3 ) = (1 + V22z% + z 3 )(l - v^zi + z 3 ). 

This means each Q-eigenvalue has its roots consisting of one square-root for both 10 + 3\/T[ 
and 10 — 3\/TT, which are zeros of the Drinfeld polynomial, as conjectured in section 3.2 of [29] 
about the relation between roots of the Q-eigenvalue and the Drinfeld polynomial in the general 
situation. Note that the above example also shows that there are eigenvectors of Q(s) which are 
not eigenvectors of the operator S = ® e erf. 
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4.3 Mathematical verification of functional relations for a general L in the root- 
of-unity six-vertex model 

In this subsection, we give a mathematical proof of the following theorem about functional relations 
of the six- vertex model at roots of unity. 

Theorem 4.2 The functional relations (3.21) and (3.22) hold for the Q-operator defined in (4.20). 
By Theorem 3.1, it suffices to show the relation (3.21) for j = N, which is equivalent to 

NQ(s)Q(-s) =T^ N \u;t). (4.31) 
By (4.20), the matrix NQ(s)Q(—s) on V is given by 
s- L (a 1 ,...,a L \NQ(s)Q(-s)\(3 1 ,...,p L ) 

= £fc= E 7< =±i</ 2 ^9 5 ? 3 (-s) 5 (-1) 5 



where ^ = £^ Pi = (mod iV). Denote 

w k (a,p-s) = (-l)^( q - k -^(-s)^- q k H-^ { - s) ^ (4.32) 
for £ Zjv, ct,P = ±1. Then we obtain the following entry-expression of the matrix NQ(s)Q(—s): 
(«!,.. . ,a L |jVQ(s)Q(-s)|/3i, ...,Pl) = s L J2k=o IlLi w fc| ^f-i (<*<,/%; s) (4.33) 

In the definition of six-vertex transfer matrix (2.2), L(s) in (2.1) is the L-operator with the C 2 - 
auxiliary space and the C 2 -quantum-space operator-entries, expressed in terms of a(s),6(s),c in 
(4.1) with the relation b(s) = a(qs). Similarly, the fusion matrix T^\t) can be obtained as the trace 
of L^\s), the fused L-operator with the C 2 -quantum space and C J -auxiliary space associated to 
L(s), where the G 7 -auxiliary space is the space of completely symmetric (j — l)-tensors of the 

C 2 -auxiliary space with the projection : <g> C 2 — > G 7 . The L (^-operator, constructed from 
the Yang-Baxter solution P^(L(s) ® a ux L(qs) <8> ■ ■ ■ ® a ux L(q J ~ 2 s)) on the C J -auxiliary space, has 
been extensively studied in the six-vertex model [21, 22, 24, 26, 30]. In order to exhibit the T^- 
operator discussed in this paper, we now employ the fusion construction to compute the explicit 
form of L^\s). Using the canonical basis {ejj^}^~Q of C J -auxiliary space associated to the C 2 -basis 
e± = | ± 1), and its dual basis {e^*Y k Z} , 

e k^ := e + 1 ke - = 77=T) £/3i+...+/3 J _ 1 =j-2fc \Pl, ■■ ■ )/^'-l)> 
e k = £ai+...+a ;( _i=j-2fc( a l) • • • i 



one can write L^\s) = (l^J(s)\ where L^j(s) is the C -(quantum-space) operator 

V ' / 0<k,l<j-l 

expressed by 



T (j) , , _ (e^*\L(s) ® aux L(qs) ® aux ■ ■ ■ ® aux L(qJ 2 s)|ep ) ) 
k,l\ s ) ~ rrJ-3" -■ ^ 



nun?*) 



21 



Consider the auxiliary-space tensor C J <8> C 2 as a subspace of (8) C 2 with the identification 



Denote := ejf } ®e_- ejjji <g> e+ for jfe = 0, . . . , j - 2. Then e[ j+1) , form a basis of 

C J ' ® C 2 with the dual basis ep +1) *, f { k j ~ 1] * given by 



e 



(j+i)* _ (j)* ^ /_| _,_ Jj)* 
fc+i ~~ e 



® (-1 + e££®(+|, 



One has 

= (4 J+1) 1^ +1) ( S )|ep +1) ) = ^y(4 J ' +1) ^ 0) ( 5 ) ®«« L (9 j_1 s)|ep +1) ). 

By induction argument, one can derive the explicit form of ]Jf\{s), and its relation with the six- 
vertex T^-operator in subsection 3.1. Furthermore the fusion relation of T^-operators follows 
from the equalities: 

(4 j+1) *\L^( S ) ®aux L{ q ^ S )\ft 1] ) = 0, 

(ft 1] *\L ij \s) ®au X L( q ^s)\ft 1] ) = b^s)(4-^\L^{s)\^). 
Therefore we obtain the following result: 

Proposition 4.2 For an integer j > 2, the C 2 -operators L^](s) (0 < k, I < j — 1) are zeros except 
k — I = 0, ±1, in which cases L^](= L k ](s)) are expressed by 



{j) _(a(q k s) \ {j) (0 - \ {j) ( 

k,k I Q a{q j-i-k s) J > h k+i,k I Q Q h ^fc-i,* I ?fc _ g _ fc Q 



The fusion matrix T^\t) is given by 

T^\t) = (s^) L t W ((g) Lf(s)), for s € C, j > 2, (4.34) 

£=1 



which satisfy the fusion relations (3.9) and (3.11) for r = 0. 
□ 

Remark. (1) Proposition 4.2, except the boundary fusion relation (3.11) valid only for the iVth 

root-of-unity q, holds as well for an arbitrary q. 

(2) The L-operator for the fusion matrix T^\t) is given by sq^2~ L^\s), whose entries are 
s-polynomials with the zero constant term (sq^ L^\(s))\ s= q except the diagonal ones: 

^ 2 -q k ~^ )' - k - j ~ 1 - 
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This implies the V-operator T (j \0) = j • Id v . □ 

We now show the Nth QQ-relation (4.31) by using the explicit form of L^ N \s) in Proposition 
4.2. By (4.34) and q N = 1, T( N \uit) is equal to the trace of the Lth monodromy matrix for the 
local operator sL^ N \qs), 

L 

T^ N \ui) = s L tr cN ((g)4 N \qs)). (4.35) 
e=i 

By Proposition 4.2, L^\qs) are zero except k — I = 0, ±1: 



4>s>=r* (1 n 1;s) ° . ), 

\ iu fc (-l,-l;s) j 

L k N k-i(l s ) = ( n ^^'n J ' L k N k+i(l s ) =( ( i , x n 

V y V io fc (-l,l;s) 



(4.36) 



where ^(a, /?; s) are functions defined in (4.32). Note that the Wk{a, (3; s) in the above expressions 
automatically determines the index (k,l) of L^(gs) where it appears, by the relation (A;, I) = 
(k,k+^). For a given pair a,j3 = ±1, the non-vanishing L^\qs) having an entry Wk(a, (3; s) 
are indexed by k £ Zjy. When evaluating the entry (a±, . . . , ctL\T( N \ujt)\f3i, . . . , (3l) in (4.35), 
the non-zero contribution of (Ll^)^ ; ^(gs) at the site I is related to that at the first site k := 
k\ £ Z^v by the relation kg = k + 2^f=i ^ i ~£ Li ■ Hence using the Wk(a, (3; s)-expression of the 
non-zero L^\qs) in (4.36), one finds (cci, . . . ,aL\T^ N \ut)\/3i, . . . with the same expression 
as (a\, ... ,aL\NQ(s)Q(— s)\(3\, (3l) in (4.33). Therefore, we obtain the relation (4.31), then 
follows Theorem 4.2. 



5 Concluding Remarks 

As a parallel theory to the Onsager-algebra symmetry of the superintegrable CPM, we have derived 
the set of functional equations of the root-of-unity six-vertex model with the Q-operator encoding 
the loop-algebra symmetry of the model. By the similar construction of the Q-operator for the 
root-of-unity eight-vertex in [3], we obtain the explicit form of Q-operator of the six- vertex model 
at roots of unity, and present the work in a way from the functional-relation aspect. As a check 
on our reasoning, we have first verified the Q-functional relation for the size L from 2 to 6 by 
direct computations, then proceed to justify mathematically that the functional relations hold for 
the six-vertex Q-operator constructed here for a general L, much as in the case of superintegrable 
CPM. These results about the Q-operator of the six-vertex model are mainly derived from the 
mathematical structures it shares commonly with the chiral Potts transfer matrix in CPM. However 
a possible physical interpretation of the root-of-vertex six-vertex Q-operator is lacking at present. 
As suggested by the role of Q-operator in CPM, the problem appears to be an interesting one, and 
the answer could be valuable as well for the symmetry problems of other solvable lattice models, 
e.g., the root-of-unity eight-vertex model in [19]. Note that in the study of the root-of-unity six- 
vertex model, the chain size L, the root's order N, and the spin S z can be arbitrary in general 
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(e.g., in [16, 17, 20]). But, just to keep things simple, we have in this paper restricted our attention 
only to even L, odd N and S z = (mod N), and leave possible generalizations to future work. 
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